April 1995 



INP Report No. 1670/PH 



QCD Analysis of Deep Inelastic Diffractive Scattering 

at HERA 

■ - - , K. Golec-Biernat and J. Kwieciriski 

in . 

Os ■ 

■ Department of Theoretical Physics, H. Niewodniczariski Institute of Nuclear Physics, 

ul. Radzikowskiego 152, 31-342 Krakow, Poland 

<• 

in : 

t— I ■ Abstract 
> 

The QCD analysis of deep inelastic diffractive scattering is performed assuming 
. the dominance of the "soft" pomeron exchange and simple, physically motivated 

parametrization of parton distributions in a pomeron. Results of the analysis are 
' compared with the recent data obtained by the HI collaboration at HERA. Both 

■ the LO and NLO approximations are considered and the theoretical predictions con- 
! cerning the quantity R = /F® for diffractive structure functions are presented. 
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The diffractive deep inelastic scattering is the process 

><■ e + p^e' + X+p' (1) 

where there is a large rapidity gap between the diffractively produced system X and 
the recoiled proton p' (or excited proton state). To be precise, the process (|1|) reflects 
the diffractive interaction of the virtual photon (see Fig. la). The recent measurements 
performed by both ZEUS and HI collaborations at the DESY ep collider HERA |3| 
have shown that the diffractive structure function constitutes a significant part of the 
proton structure function F 2 . 

The cross section for the process (1) is related in the following way to the diffractive 
structure functions dF ^f,t) <(; f yt) ; 

dxpdt dxpdt 
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where the kinematical variables are defined below: 

PQ 

PeP 

-o2 2 Q 

Q = - q x = — 
2pq 

t = (p-p') 2 /3= 9 ( Q2 , x P = X - (3) 
2q(p - pf) f3 

and where q = p e — p' e . The four momenta p e ,p' e ,P and p' correspond to those of the ini- 
tial electron, final electron, initial proton and the recoiled proton respectively (see Fig. 1). 

The main aim of this paper is to analyse the diffractive processes in deep inelastic 
ep scattering assuming that they are dominated by a "soft" pomeron exchange J|- ||10|| . 



Different description of this processes based on the "hard" pomeron which follows from 
perturbative QCD has been discussed in [11] -[14]. The "hard" pomeron would however 



lead to much steeper xp dependence than that which is observed experimentally ||. 

Assuming that the diffractive process , y*{Q 2 )+p — > X+p is dominated by the pomeron 
exchange (see Fig. lb) with the pomeron being described as a Regge pole with its trajec- 
tory 

a P (t) = ctp(O) +a' P t (4) 

where for the "soft" pomeron «p(0) ~ 0.08, a' P = 0.25 GeV -2 , we get the following 
factorizable expression for the diffractive structure functions : 

dF 2 D L (x P ,(3,Q\t) ti ^ ^ p 



dxpdt 



f(x P ,t) F*(/3,Q\t) (5) 



The function f(xp,t) is the so called "pomeron flux factor" which, if the diffractively 
recoiled system is a single proton, has the following form |10| : 



f{x P ,t) =N x P —— (6) 

lD7T 

where B(t) describes the pomeron coupling to a proton. Its normalization is such that 
the pomeron contribution to the pp total cross-section <J pp (s) is 

/ s N Q P (0)-1 



a„{8) = B i QS) - (7) 

with s = lGeV 2 . The normalization factor iV was set to be equal to - following the 
convention of refs. 0, |6|, FToT . We have also set UTOl 



B(t) = 4.6 mb 1/2 exp(1.9GeV~ 2 t) 

and have slightly increased a P (0) in the flux factor f(xp,t) (J[^) up to «p(0) = 1.1 
The slope of the pomeron trajectory was kept equal to a' P = 0.25 GeV -2 . 
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The functions F[ L (f3,Q 2 ,t) are the pomeron structure functions with the variable (3 
playing the role of the Bjorken scaling variable for the 7*(Q 2 ) pomeron inelastic "scat- 
tering". The factorization property of the diffractive structure function into the pomeron 
"flux" and the pomeron structure function is the direct consequence of the single pomeron 
exchange mechanism of diffraction and of the assumption that the pomeron is described 
by a single Regge pole. Factorization does not in general hold in models which are based 
entirely on perturbative QCD [T(], |I5| . Recent measurements || have however found that 
the factorization of the diffractive structure functions agrees very well with the experi- 
mental data. It has moreover been established that the xp dependence of the flux factor 
(integrated over t) can be parametrized as Xp n with n = 1.19. This observation confirms 
dominance of the soft pomeron exchange since the "hard" pomeron with its intercept 
op(0) ~ 1.5 would give much steeper behaviour ~ Xp 2 or so. 



In the region of large Q 2 the pomeron structure function F 2 P is expected to be described 
by the "hand-bag" diagram (see Fig. 2a) [|], ||. It provides the partonic description of deep 
inelastic diffraction which leads to the Bjorken scaling mildly violated by the logarithmic 
scaling violations implied by perturbative QCD. In the QCD improved parton model 
the pomeron structure function F 2 p (/3, Q 2 , t) is related in conventional way to the quark 
qf(/3,Q 2 ,t) distributions in a pomeron: 

Ff(P,Q 2 ,t) = 2p £e 2 gf(/?,Q 2 ,t) (9) 

i 

where ti are the quark charges (note that q[(/3,Q 2 ,t) = qf(/3,Q 2 ,t)). The formula (0) 
acquires 0(a s ) corrections in the next-to-leading approximation of perturbative QCD. 
The next-to-leading approximation is also needed for theoretically consistent introduc- 
tion of the pomeron longitudinal structure function F[(f3, Q 2 ,t). In what follows we shall 
therefore discuss both the leading and next to leading approximations. In the leading 
order approximation we shall set F[ = 0. 



At first we shall specify the details of the parton distributions in the pomeron at the 
reference scale Q 2 = Ql which we set equal to 4GeV 2 . 

At small (3 both the quark and gluon distributions are assumed to be dominated by 
the pomeron exchange as illustrated in the Fig. 2b. 

Pg p (P,Ql,t) = a g (t)^- a ^ (10) 

The functions a«(t) and a g (t) can be estimated from the factorization of pomeron couplings 
1J,1 (see Fig. 2b): 

di(t) = r(t) af 

a g (t) = r(t) d> g (11) 
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where the parameters af and a p g are the pomeron couplings controlling the normalization 
of the small x behaviour of the sea quark and gluon distributions in a proton i.e. 



xqf(x, Qo) + xqf(x } Q 



i 



x<f(x,Ql) = a?x l - a ^ (12) 



and the function r(t) is: 



= ^Gpppjt) 

v 7 2 5(0) V ; 

The coupling Gppp(t) is the triple pomeron coupling (see Fig. 2 b) and its magnitude can 
be estimated from the cross-section of the diffractive production p + p — > p + X in the 
limit of large mass Mx of the diffractively produced system X. The relevant formula for 
this cross-section is: 

(The factor | in the eq. ( P~3| ) follows from the convention iV = - for the normalization 
factor N in the eq. @ 0, || |10|1). We neglected the (weak) t dependence of the function 
r(£) and have estimated its magintude from the recent Tevatron data [IB as r(t) ~ r(0) = 



0.089. The parameters af were estimated assuming that the sea quark distributions in a 
proton can be parametrized as: 

xqf(x, Ql) + xqf(x, Q 2 ) = 2 af x 1 ""^ ) (1 - x) 7 (15) 

and fixing the constants a\ from the requirement that the average momentum fraction 
which corresponds to those distributions is the same as that which follows from the recent 



parametrization of parton distributions in a proton [|T7| , [18). The momentum sum rule 
has also been used to fix the parameter a v g i.e. we assumed 

xg p {x, Ql) = aP g x 1 ' ^ (1 - xf (16) 

and imposed the condition that the gluons carry 1/2 momentum of the proton. We 
extrapolated the pomeron dominated quark and gluon distributions in a pomeron (see 
flT0|) ) to the region of arbitrary values of (3 by multiplying the factor p 1 ~ a p(°) by 1 — (3 0. 

We have also included the term proportional to /3(1 — /3) in both the quark and gluon 
distributions 0. The normalization of this term in the quark distributions has been 
estimated in || assuming that it is dominated by the quark-box diagram of Fig. 2c with 
the non-perturbative couplings of pomeron to quarks. In this model one gets: 

Pq p (P,Qlt) = ^P(l-P) (17) 

where C ~ 0.17 |J. We found that the fairly reasonable description of data can be 
achieved provided that the constant C is enhanced by a factor equal to 1.5. We have 
also assumed that the relative normalization of the quark distributions in a pomeron 
corresponding to different flavours is the same as that of the sea quark distributions in 
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a proton [17, 18[j. Finally the normalization of the term proportional to f3(l — /3) in the 
gluon distribution in a pomeron has been obtained by imposing the momentum sum rule. 

As the result of the estimates and extrapolations discussed above the parametrization 
of parton distributions in a pomeron at the reference scale Ql = 4GeV 2 looks as follows: 



Pu p (P,Ql,t) 
(3 d p ((3,Q 2 ,t) 
(3s p (f3,Qlt) 
(3 c p ((3,Q 2 ,t) 
(3 g p ((3,Q 2 ,t) 



;i - s) s p ((3) 
;i - s) s p (p) 



0.4 

0.4 > 
0.2 (1 - 5) S p {f3) 
5 S p {(3) 

(0.218 /r a08 + 3.30 0) (1-/3) 



where the function S p (f3) is parametrized as below 



S F (f3) = (0.0528 p~ om + 0.801 (3) (1 - (3) 



(19) 



and 5=0.02 jlT], Following the approximations discussed above we have neglected 
the t dependence in those parton distributions. The analysis of the pomeron structure 
functions based on different parametrizations of parton distributions in a pomeron has 
recently been presented in refs. || [L9| . 

The parton distributions defined by eqs.(0,0) were next evolved up to the values of 



Q for which the data exist using the LO Altarelli-Parisi evolution equations [20, |2T[ with 



A = 0.255 GeV [|11|. In Fig. 3a,b we show our results for the quantity: 



dxi 



XPL 



dt dF 2 D (x P ,P,Q 2 ,t) 
dx pdt 



(20) 



with xpl = 0.0003 and xph = 0.05 0. The theoretical curves are compared with the 
recent data from the HI collaboration at HERA H. We find that the agreement with the 
data is very good concerning both the shape in (3 as well as the evolution with Q 2 . The 
following points should be noticed: 



• The diffractive structure function F^(/3, Q 2 )at the moderate values of (3 is dom- 
inated by "hard" component of the quark distributions in a pomeron (i.e. that 
component which is proportional to (3(1 — (3)). Relatively small magnitude of the 
soft component (i.e. of the term which is proportional to (3~ om ) is a direct conse- 
quence of the small magnitude of the triple pomeron coupling (see the eqs. (|T0"D , (|Tl"D , 
(0) and ©). 

• The soft component of the parton distributions is important at the small values of 
(3 < 0.1 or so. The agreement of our parametrization with the data shows that the 
estimate of its normalization from the factorization of pomeron couplings is quite 
reasonable. At large Q 2 the f3 spectrum becomes softer as the result of the QCD 
evolution. 
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• It turns out that gluons carry more than 60 % of the pomeron momentum. The 
gluon content of the pomeron can only indirectly manifest itself through the QCD 
evolution of the diffractive structure function with Q 2 . It should be observed at this 
point that the relatively large gluon distributions in a pomeron with a "hard" 1 — /3 
spectrum leads to the increase of the diffractive structure functions with increasing 
Q 2 for fixed (3 , up to rather relatively large values of /3 ~ 0.5. It may be seen from 
Fig. 3a that this effect is nicely confirmed by the data. It should be noted at this 
point that the proton structure function F%(x, Q 2 ) starts to decrease with increasing 
Q 2 already for x > 0.1. 



We have also performed the NLO QCD analysis in order to be able to consistently 
introduce the quantity R = F^(f3,Q 2 )/Ff(f3 } Q 2 ) where F^(f3,Q 2 ) = F°(f3,Q 2 ) - 
F°(f3,Q 2 ) and 

^(ftcn-rW (21) 

Jxpl J— oo dxpdt 

We compare the NLO and LO results for F®(/3, Q 2 ) in Fig. 4 a,b. We can see that both 
approximations lead to similar predictions. In Fig. 4 c,d we show our results for R. The 
longitudinal diffractive structure function is mainly driven by the gluon distributions in a 
pomeron and a large amount of gluons (see [18]) implies that R can reach 0.5 for (3 < 0.1. 
It may also be seen from Figs. 4 a,b that R exhibits much more softer dependence on (3 
than the structure function F 2 p (/3, Q 2 ). 



To summarize, we have presented in this paper the QCD analysis of the diffractive 
structure function within the simple soft pomeron exchange picture. In this model the 
diffractive structure functions factorize into the flux factor and the pomeron structure 
functions which are sensitive on the parton content of the pomeron. We based our anal- 
ysis on the simple parametrization of the parton distributions which utilized both the 
factorization of pomeron couplings and the quark-box contribution and obtained a good 
description of the experimental data. We have also performed the NLO QCD analysis 
and estimated the magnitude of R = F® (/3, Q 2 )/F^(/3, Q 2 ) (see Fig. 4 c,d). Its knowledge 
should be useful for more refined experimental determination of the diffractive structure 
function F 2 D . 
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Figure Captions 



1. (a) Kinematics of the large rapidity gap process e + p^e' + X+ p'.The wavy line 
represents the virtual photon. 

(b)The Pomeron exchange diagram for diffractive production of the hadronic system 
X by a virtual photon. The wavy and zigzag lines represent the virtual photon and 
pomeron respectively. 

2. (a)The hand bag diagram for the virtual Compton diffractive production. The 
wavy and zigzag lines represent the virtual photon and pomeron respectively. The 
continuous lines in the upper part of the diagram denote quarks (antiquarks). 

(b) The triple pomeron diagram contribution to the "hand-bag" diagram of Fig. 2a. 
The coupling Gppp denotes the triple pomeron coupling. 

(c) The quark box diagram contribution to the "hand-bag" diagram of Fig. 2a. 

3. Theoretical predictions for the diffractive structure function F 2 D (/3,Q 2 ) defined by 
the formula ( ^0|) and their comparison with the data from HERA ||. The structure 
function -F 2 D (/5, Q 2 ) is plotted (a) as the function of Q 2 for fixed values of (3 and (b) 
as the function of (3 for fixed values of Q 2 . 

4. (a) Comparison of LO (continuous lines) with NLO (dashed lines) results for the 
structure function F 2 D (/5,Q 2 ) defined by the eq. (|20|) plotted as the function of Q 2 
for fixed values of (3. 

(b) Comparison of LO (continuous lines) with NLO (dashed lines) results for the 
structure function F®(j3,Q 2 ) defined by the eq. (|20|) plotted as the function of (3 
for fixed values of Q 2 . 



5. (a) Theoretical predictions for the quantity R 
for fixed values of (3. 

(b) Theoretical predictions for the quantity R ■ 
fixed values of Q 2 . 




plotted as the function of Q 2 




plotted as the function of (3 for 
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Figure 3 




Figure 4 




Figure 5 




